Superfluidity in many fermion systems: Exact renormalisation group treatment 
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The application of the exact renormalisation group to symmetric as well as asymmetric many- 
fermion systems with a short-range attractive force is studied. Assuming an ansatz for the effective 
action with effective bosons, describing pairing effects and a set of approximate flow equations for the 
effective coupling including boson and fermionic fluctuations has been derived. The phase transition 
to a phase with broken symmetry is found at a critical value of the running scale. The mean-field 
results are recovered if boson-loop effects are omitted. The calculations with two different forms 
of the regulator are shown to lead to a similar results. We find that, being quite small in the case 
of the symmetric many-fermion system the corrections to mean field approximation becomes more 
important with increasing mass asymmetry. 



I. INTRODUCTION 

There is a growing interest in applying the Exact 
Renormalisation Group (ERG) formalism to few- and 
many-body systems 0, 0, Q when the underlying inter- 
action is essentially non-perturbative. Regardless of the 
details all ERG-based approaches share the same distinc- 
tive feature, a successive elimination/suppression of some 
modes, resulting in effective interaction between the re- 
maining degrees of freedom. One specific way of imple- 
menting such a procedure is to eliminate modes by ap- 
plying a momentum-space blocking transformation with 
some physically motivated cutoff. The effect of varying 
a cutoff is described by nonlinear ERG evolution equa- 
tions, which include the effect of the eliminated modes. In 
the following we will use the variant of ERG based on 
the concept of Average Effective Action (AEA) 4] . The 
corresponding evolution equation can written in the fol- 
lowing general form 



Tr 



(1) 



Here I^ 2 ) is the second functional derivative of the AEA 
taken with respect to all types of field included in the 
action, and R is a regulator which should suppress the 
contributions of states with momenta less than or of the 
order of running scale k. To recover the full effective ac- 
tion we require R(k) to vanish as fc — > 0, in other respects 
its form is rather arbitrary. The concrete functional form 
of the regulator has no effect on physical results pro- 
vided no approximations/truncations were made. By solv- 
ing the ERG equations one can find a scale dependence 
of the coupling constants and thus determine the ERG 
flow. The whole approach is nonperturbative so that 
some physically motivated assumption about the func- 
tional form of the effective action should be made It is 
well known that the many-fermion system with attrac- 
tive interaction favours the formation of the correlated 
fermion pairs leading to the symmetry breaking and re- 
lated phenomena. Since we expect the appearance of the 
correlated fermion pairs in a physical ground state, we 
need to parametrise our effective action in a way that 



can describe the qualitative change in the physics when 
this occurs. A natural way to do this is to introduce a 
boson field whose vacuum expectation value (VEV) de- 
scribes this correlated pair Q and study the evolution of 
this effective degrees of freedom. At the start of the RG 
evolution, the boson field is not dynamical and is intro- 
duced through a Hubbard-Stratonovich transformation 
of the four-point interaction. As we integrate out more 
and more of the fermion degrees of freedom by running 
the cutoff scale k to lower values, we generate dynami- 
cal terms in the bosonic effective action. In this paper 
we treat both symmetric and asymmetric many fermion 
systems. The corresponding ansatz for the boson-fermion 
effective action consists of the kinetic terms for boson and 
fcrmions and the interaction term and for two types of 
fermions can be written as 



T\fi,k] 



d 4 x(r B [ M ,fe] + rir[/i,fc] + rj[fc]) (2) 



Here ^b(f) is the boson (fermion) part of AEA 



^ (Zf (iOt +Va+Vb) + ^ V 2 ) Cj> - U(<j>, f ) 

(3) 



r> = ]T 4 [z^idt + Mi ) + ^ v 2 ) <fc (4) 

i—a ^ 

and Ti is the interaction term 



(5) 



M is the reduced mass of the fermion in vacuum and 
the factor l/2m with m = M a + M& in the boson ki- 
netic term is chosen simply to make Z m dimensionless. 
The coupling Z g , the wave-function renormalisations fac- 
tors Zrf,^ and the kinetic-mass renormalisations factors 
Z m ,M all run with on k, the scale of the regulator. Hav- 
ing in mind the future applications to the crossover from 
BCS to BEC (where chemical potential becomes nega- 
tive) we also let the chemical potentials [i a and fib run, 



f 



thus keeping the corresponding densities (and Fermi mo- 
menta PF,i) constant. The bosons are , in principle, cou- 
pled to the chemical potentials via a quadratic term in <j>, 
but this can be absorbed into the potential by defining 
U = U — (pi + p<2)Z$$ <f>. The evolution equations in- 
clude running of chemical potentials, effective potential 
and all couplings (Z^, Z m , Zm,i, Z^^, Z g ). However, in 
this paper we allow to run only Z^, parameters in the 
effective potential (u's and po) and chemical potentials 
since this is the minimal set needed to include the ef- 
fective boson dynamics. The system with one type of 
fcrmion corresponds to the limit M a = Mb. 

We expand the effective potential about its minimum, 
(jy(j) = po,so that the coefficients itj are defined at p = po, 

U(p) = u +u 1 (p-p ) + ^ u 2 {p-pof + ^ u 3 (p-p Q ) 3 + - ■ ■ , 

(6) 

where we have introduced p = <fy<f>. Similar expansion 
can be written for the renormalisation factors. The co- 
efficients of the expansion run with the scale. The phase 
of the system is determined by the coefficient u\. We 
start evolution at high scale where the system is in the 
symmetric phase so that u\ > 0. When the running 
scale becomes comparable with the pairing scale (close 
to average Fermi-momentum) the system undergoes the 
phase transition to the phase with broken symmetry, en- 
ergy gap etc. The point of the transition corresponds to 
the scale where u\ = 0. The bosonic excitations in the 
gapped phase are gap-less Goldstone bosons. Note, that 
in this phase the minimum of the potential will also run 
with the scale k so that the value po(k — > 0) determines 
the physical gap. 

The important part of any ERG treatment is the choice 
of the regulator. Ideally, the physical results should not 
depend on this choice. However, some sort of trunca- 
tions and approximations should always be made in real 
calculations to render the system of the resulting evolu- 
tion equations solvable so that the convenient choice of 
the regulator is the question of significant practical im- 
portance. In our approach the boson regulator has the 
structure 



Rb = Rsdiag(l, 1), 



(7) 



and the fermion regulator for both types of fermions has 
the structure 

R F ,i = sgn(ei(q) - p^Rp.^q, pi,k)diag(l, -1) (8) 

Note that this regulator is positive for particle states 
above the Fermi surface and negative for the hole states 
below the Fermi surface. The function Rp should sup- 
press the contributions of states with momenta near the 
Fermi surface, \q — pf\ ~ k. Once a a large gap has ap- 
peared in the fermion spectrum, there are no low-energy 
fermion excitations and so the fermionic regulator plays 
little further role. However, while the gap is zero or small, 
it is crucial that the sign of the regulator match that of 



the energy, q 2 /2M — p, and hence it is p, which appears 
in the sign functions. 

The other important part of the ERG approach is fix- 
ing the boundary conditions which define the form of the 
AEA at some initial scale so that at some large starting 
scale k = K we demand that the Lagrangian be equiva- 
lent to a purely fermionic theory with the contact inter- 
action 



(9) 



Here Cq(K) is the strength of the energy- independent 
term in the effective NN interaction in vacuum. This is 
evaluated at the scale K, using the same regularisation 
procedure as we apply in matter. This equation implies 
that u\ and g at the this scale are related by 



C (K) 



ui{K)' 



(10) 



Using this boundary condition we can relate the pair- 
ing phenomena in a many body invironment with the 
fermion-fermion interaction in vacuum. 

The boson potential U is obtained by evaluating the 
effective action for uniform boson fields. It evolves ac- 
cording to 



d k U 



1 

Va 



d k T, 



(11) 



where V4 is the volume of spacetime. Substituting our 
expansion of U, Eq. JJjJl, on the left-hand side leads to a 
set of ordinary differential equations for the u n . 

For the boson wave- function renormalisation factor, 
Z$, we need to consider a time-dependent background 
field. Taking 



-ipot 



(12) 



where r\ is a constant, we can get the evolution of Z^ 
from 



z 1 d ( 52 d r 

k V4 dp ydrjdrf k 



7J = 



(13) 



Po=0 



Let us first consider the results of the calculations in the 
case of symmetric fermion matter. We solve the evolution 
equations numerically with two types of cutoffs. First, we 
use the smoothed step-function type of regulator (called 
hereafter as Rif) 

k 2 k 2 

RlF = 7777^1(9 ~ Pf, k, a); Rib = 77—^1(9, k, a) (14) 
ZM Zm 



where 

6 1 (q,k,a) 



2er/(l/er) 



erf{^— ) + erf{— — ) 



ka 



ka 



(15) 
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FIG. 1: Numerical solutions to the evolution equations for 
infinite ao and pf = 1.37 fm, starting from K — 16 frn -1 . We 
show the evolution of all relevant parameters for the cases of 
fermion loops only (orange/grey lines), and of bosonic loops 
with a running (blue/black lines). All quantities are ex- 
pressed in appropriate powers of fm^ 1 . 

with a being a parameter determining the sharpness of 
the step. Second, we use the sharp cutoff function i?2_F 
which is somewhat similar to one suggested in Rcf. 6] for 
the pure boson case and chosen in a rather peculiar way 
to make the calculations as simple as possible 

R2F = 7^i(( k +P») 2 - <1 2 )0(P, + k-q) + 

((* + prf + q 2 - 2pl)9(q - Pfi + k)], (16) 

R2B = -^(k 2 -q 2 )6(k-q), (17) 

where p^ = (2M/^) 1 / 2 . The fermion sharp cutoff consists 
of two terms which result in modification of the parti- 
cle and hole propagators respectively. The hole term is 
further modified to suppress the contribution from the 
surface terms, which may bring in the dangerous depen- 
dence of the regulator on the cutoff scale even at the van- 
ishingly small k. As an example, we focus on the param- 
eters relevant to neutron matter: M — 4.76 fm~ 1 ,pp — 
1.37/m . Let's first discuss the results obtained with 
the smooth cutoff R\. We found that the value of the 
physical gap is practically independent of either the val- 
ues of the width parameter a (varied within some range) 
or the starting scale K provided K > 5/m _1 . The re- 
sults of the calculations are shown on Fig. 1 At starting 
scale the system is in the symmetric phase and remains 
in this phase until u\ hits zero at k cr u ~ 1.2/m _1 where 
the artificial second order phase transition to a broken 
phase occurs and the energy gap is formed. Already at 
k ~ 0.5 the running scale has essentually no effect on the 
gap. It is worth mentioning that we found very small (on 
the level of 1%) contribution to the gap from the boson 
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FIG. 2: Evolution of the parameters when the sharp cut-off 
is used 



loops, due to cancellations between the direct contribu- 
tions to the running of the gap and indirect ones via u^. 
The boson loops play much more important role in the 
evolution of it 2 and Z$. In fact, they drive both couplings 
to zero at k — > 0. We note however, that the effect of the 
boson loops for the gap may still be more visible if the 
evolution of the other couplings is included. 
The results obtained with sharp cutoff regulator which 
are shown on Fig. 2. One immediate observation is that 
that the results become starting scale independent as 
long as K > 5/m _1 similar to the case with the smooth 
cutoff. However, the artificial phase transition occurs at 
lower values of the running scale k ~ 0.7 fm. At ap- 
proximately k ~ 0.2/to _1 the value of the gap becomes 
scale independent. One notes that the curves obtained 
with different regulators and describing the evolution of 
the gap, being rather different at intermediate scales, ap- 
proach each other with decreasing scale resulting in very 
close values for the physical gap. This is an encouraging 
results taking into account that, although the hypotheti- 
cal exact results must be independent of the choice of the 
regulator in practice it is not garanteed at all given the 
assumed ansatz for the effective action and truncations 
made. The same conclusion also holds for the other quan- 
tities. The couplings Z$ and 1*2 first grow with scale and 
then start decreasing eventually coming to zero. Chemi- 
cal potential begins to decrease at the point of phase tran- 
sition and becomes scale independent at k ~ 0.2/m _1 . 
However, in this case the numerical values of the chemi- 
cal potentials obtained with different regulators differ by 
approximately 20% so that this quantity is more sensitive 
to the details of effective action and to the trancations 
made. 

Let us now discuss the results for the many fermion 
system with two fermion species. For simplicity we con- 
sider the case of the hypothetical "nuclear" matter with 
short range attractive interaction between two types of 
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FIG. 3: Evolution of the gap in the MF approach (dashed 
curve) and with boson loops (solid curve) in the unitary 
regime a = — oo as a function of a mass asymmetry. 



fcrmions, light and heavy, and study the behaviour of 
the energy gap as the function of the mass asymmetry. 
We choose the Fermi momentum to be pp = 1.37/to -1 . 
One notes that the formalism is applicable to any type 
of a many-body system with two fermion species from 
quark matter to fermionic atoms so that the hypothet- 
ical asymmetrical "nuclear" matter is simply chosen as 
a study case. We assume that M a < Mb, where M a is 
always the mass of the physical nucleon. 

First we consider the case of the unitary limit with the 
infinite scattering length. The results of our calculations 
for the gap are shown on Fig. 3. 

We see from this figure that increasing mass asymme- 
try leads to a decreasing gap that seems to be a natural 
result. However, the effect of the boson loops is found 
to be small. We found essentially no effect in symmet- 
ric phase, 2 — 4% corrections for the value of the gap 
in the broken phase and even smaller corrections for the 
chemical potential so that one can conclude that the MF 
approach indeed provides the reliable description in the 
unitary limit for both small and large mass asymmetries. 
It is worth mentioning that, similar to above considered 
case of the fermion matter with one type of fermion, the 
boson contributions are more important for the evolu- 
tion of U2 where they drive u 2 to zero as k — > making 
the effective potential convex in agreement with the gen- 
eral expectations. This tendency retains in the unitary 
regime regardless of the mass asymmetry. 

We have also considered the behaviour of the gap as 
the function of the parameter ppa for the cases of the 
zero asymmetry M a — Mb and the maximal asymmetry 
Mb = WM a . The results are shown on Fig. 4. 

One can see from Fig. 4 that in the case of zero (or 
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FIG. 4: Evolution of the gap as a function of the parameter 
PFa. The upper pair of the curves corresponds to the calcula- 
tions with no asymmetry in the MF approach (dashed curve) 
and with boson loops (solid curve) and the lower pair of the 
curves describes the results of calculations with the maximal 
asymmetry when Mt = 10M a 



small) asymmetry the corrections stemming from boson 
loops are small at all values of the parameter ppa con- 
sidered here (down to ppa = 0.94). On the contrary, 
when Mb = 10M a these corrections, being rather small 
at ppa > 2 becomes significant (~ 40%) when the value 
of pi? a decreases down to ppa ~ 1. We found that at 
ppa ~ 1 the effect of boson fluctuations becomes <~ 10% 
already for Mb = 5M a . One can therefore conclude that 
the regime of large mass asymmetries, which starts ap- 
proximately at Mb > 5M a , moderate scattering length 
and/or the Fermi momenta is the one where the MF de- 
scription becomes less accurate so that the calculations 
going beyond the MFA are needed. One might expect 
that the deviation from the mean field results could even 
be stronger in a general case of a large mass asymmetry 
and the mismatched Fermi surfaces but the detailed con- 
clusion can only be drawn after the actual calculations 
are performed. 

We were not able to follow the evolution of the system 
at small gap (or small ppa) because of the non-analyticity 
of the effective action in this case. This non-analyticity 
of the effective action can explicitly be demonstrated in 
the mean-field approximation. The flow equations can 
be solved analytically in this case and one can see from 
the solution, which has a closed-form expression in terms 
of an associated Legendre function, at k = 0, that 

the fermion loops contain a term <fi^ tfilog^ <f>). It remains 
to be seen whether, within the given ansatz, the full solu- 
tion of the system of the partial differential equations for 
the effective potential and running couplings is required 
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to trace the evolution of the system in the case of small 
gaps. 

In summary, we have studied the pairing effect for the 
asymmetric fermion matter with two fermion species as 
a function of fermion mass asymmetry. We found that 
regardless of the size of the fermion mass asymmetry the 
boson loop corrections are small at large enough values 
of pfcl so that the MFA provides a consistent description 
of the pairing effect in this case. However, when ppa ~ 1 
these corrections become significant at large asymmetries 
{Mb > 5M a ) making the MFA inadequate. In this case 
it seems to be necessary to go beyond the mean field 
description. 

There are several ways where this approach can fur- 
ther be developed. In the case of asymmetric systems 
the next natural step would be to consider the case of 
the mismatched Fermi surfaces taking into account the 
possibility of formation of Sarma , mixed H,|| and/or 
LOFF phases, exploring the importance of the bo- 



son loop for the stability of those phases and applying 
the approach to the real physical systems, for example 
fermionic atoms. Work in this direction is in progress. 
The other important extension of this approach would 
be to take into account running of all couplings of the 
average effective action, use different type of cut-off func- 
tion, preferably the smooth one and inlude both particle- 
hole channel and long range forces. The three body force 
effects pj|. when the correlated pair interact with the 
unpaired fermion may also be important, especially for 
non-dilute systems. 
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